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Abstract — There  has  been  a  growing  interest  in  passive  radar 
systems  in  the  research  community  over  the  last  decade  because 
of  the  several  merits  they  offer,  including  ease  of  deployment, 
low  cost,  and  non-detectability  of  the  receivers.  During  the  same 
period,  the  idea  of  distributed  MIMO  radar  and  its  advantages 
under  the  coherent  and  non-coherent  operating  scenarios  has 
been  extensively  studied.  Keeping  these  benefits  it  mind,  in  this 
paper,  we  consider  a  UMTS-based  passive  multistatic  radar  with 
distributed  antennas.  We  compute  the  ambiguity  profiles  of  this 
radar  system  under  both  the  coherent  and  non-coherent  modes. 
The  non-coherent  processing  mode  improves  the  target  detection 
performance  by  obtaining  spatially  diverse  looks  of  the  target. 
On  the  other  hand,  coherent  processing  enhances  the  resolution 
of  target  localization.  We  use  numerical  examples  to  demonstrate 
our  analytical  results. 

Index  Terms — Ambiguity  function,  coherent  processing,  dis¬ 
tributed,  multistatic,  passive  radar,  resolution,  UMTS  signals. 

1.  Introduction 

The  concept  of  passive  radar  systems  or  passive  coherent 
location  systems  has  been  in  existence  since  the  1930s 
[l]-[3].  This  topic  has  drawn  much  research  attention  only 
during  the  last  decade  [4]-[6].  When  compared  with  eonven- 
tional  monostatie  systems  that  have  the  transmitter  and  reeeiver 
eoloeated,  these  systems  often  rely  on  signals  of  opportunity 
as  the  transmitted  waveforms  and  the  reeeivers  are  eompletely 
passive  and  loeated  away  from  the  transmitters  [6],  [7].  There¬ 
fore,  the  reeeivers  are  not  deteetable  and  the  system  to  be 
deployed  does  not  require  expensive  transmitting  hardware. 
Passive  radar  systems  do  not  have  eontrol  over  the  transmitted 
waveforms  but  reeeive  a  eopy  of  them  direetly  from  the  trans¬ 
mitter  and  this  will  be  used  as  a  referenee  signal.  Signals  of 
opportunity  ean  be  from  several  sourees  ineluding  television 
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[l]-[9]  and  audio  broadcast  signals  and  FM  radio  [10],  [11], 
satellite-based  [12],  and  mobile  communications  systems  [13]. 

In  this  paper,  we  consider  a  passive  radar  system  based  on 
universal  mobile  telecommunications  systems  (UMTS).  These 
UMTS  signals  are  known  to  provide  high  resolution  and  low 
sidelobes  and  hence  are  of  partieular  interest  to  the  passive  radar 
eommunity.  Further,  with  the  rapid  inerease  in  their  usage  glob¬ 
ally,  the  eoverage  areas  for  mobile  eommunieations  systems  has 
inereased  several-fold  in  the  reeent  years.  Earlier,  performanee 
of  mobile  eommunieations-based  passive  eoherent  loeation  sys¬ 
tems  has  been  studied  using  experiments  [13]  and  simulations 
[14].  More  reeently,  in  [15],  the  authors  derived  analytieal  ex¬ 
pressions  for  the  monostatie  and  bistatie  ambiguity  funetions 
using  these  UMTS  downlink  signals.  Further,  they  eomputed  the 
modified  Cramer-Rao  lower  bound  (CRB)  for  these  systems. 

Conventional  passive  radar  systems  operate  in  a  bistatie  eon- 
figuration.  However,  very  reeently,  there  has  been  a  growing 
interest  in  using  passive  radar  systems  in  multistatie  eonfigu- 
rations.  This  is  motivated  from  the  several  advantages  of  dis¬ 
tributed  multiple  input  multiple  output  (MIMO)  radar  systems 
demonstrated  in  the  last  deeade  [16]-[20].  When  operating  in 
a  non-eoherent  eonfiguration  [21]  these  systems  exploit  spa¬ 
tially  diverse  looks  of  the  target  to  improve  the  deteetion  per¬ 
formanee  [17].  This  goal  is  aeeomplished  by  employing  widely 
spaeed  transmit  and  reeeive  antennas.  When  the  target  reflee- 
tions  are  weak  between  partieular  transmitter-reeeiver  pair,  it 
is  highly  likely  that  these  will  be  eompensated  by  some  of  the 
other  pairs  that  have  stronger  returns.  Further,  when  there  is  per- 
feet  phase  synehronization  between  all  the  transmitters  and  re¬ 
eeivers,  MIMO  radar  ean  operate  in  a  eoherent  proeessing  mode 
to  enhanee  the  target  loealization  resolution  [17],  [22].  How¬ 
ever,  note  that  when  the  antennas  are  too  far  apart,  phase  syn¬ 
ehronization  may  not  be  aehievable  and  in  sueh  seenarios,  we 
have  to  resort  to  non-eoherent  proeessing.  It  is  natural  for  pas¬ 
sive  radar  systems  to  benefit  from  the  same  reasons  by  oper¬ 
ating  in  multistatie  eonfiguration.  In  [23],  the  authors  presented 
a  single-stage  non-eoherent  proeessing  approaeh  that  performs 
target  deteetion,  loealization,  and  de-ghosting  for  passive  mul¬ 
tistatie  radar. 

In  this  work,  we  eompute  the  ambiguity  funetions  for  passive 
multistatie  radar  systems  using  UMTS  signals  (See  also  [24]). 
The  ambiguity  funetion  is  an  important  performanee  metrie  for 
radar  systems  beeause  it  deseribes  the  global  performanee  (loeal 
aeeuraey  and  sidelobes).  We  eonsider  both  the  eoherent  and 
non-eoherent  operating  seenarios  in  the  Cartesian  domain.  The 
delays  and  Doppler  shifts  assoeiated  with  a  multistatie  system 
are  non-linear  functions  of  the  geometry.  Therefore,  the  ambi¬ 
guity  funetions  are  also  geometry-dependent.  We  use  numer- 
ieal  examples  to  study  the  aehievable  performanee  of  these  sys- 
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Fig.  I.  Passive  multistatic  radar  system  (graphic  for  antennas  from  [25]). 


terns.  The  rest  of  the  paper  is  organized  as  follows.  In  Seetion  II, 
we  present  the  transmitted  signal  model  for  passive  multistatie 
radar.  In  Seetion  III,  we  present  the  non-eoherent  measurement 
model,  followed  by  the  deteetor  and  the  ambiguity  funetion. 
Similarly,  in  Seetion  IV,  we  eonsider  the  eoherent  seenario.  In 
Seetion  V,  we  present  numerieal  examples  of  our  analytieal  re¬ 
sults.  Finally,  we  provide  eoneluding  remarks  in  Seetion  VI. 

II.  Signal  Model 

First,  we  deseribe  the  radar  setup  in  this  seetion,  followed 
by  the  target  models  under  both  the  non-eoherent  and  eoherent 
operating  modes  in  the  following  seetions.  We  eonsider  a  pas¬ 
sive  multistatie  radar  system  eonsisting  of  Mt  transmitters  and 
Mr  reeeivers  that  are  widely  separated  (see  Fig.  1).  We  assume 
that  there  is  always  a  line-of-sight  eomponent  between  all  trans- 
mitter-reeeiver  pairs.  In  Fig.  1,  we  have  also  marked  the  illumi¬ 
nated  seene  of  interest,  the  referenee  paths  and  the  target  eeho 
paths  for  a  multistatie  system  with  3  transmitters  and  2  reeeivers. 

Let  the  zth  transmitter  be  loeated  at  ti  =  [txi^tyi^tzi].  Simi¬ 
larly,  the  jth  reeeiver  is  loeated  at  We  as¬ 

sume  the  transmitters  and  reeeivers  are  stationary  but  the  results 
ean  be  easily  extended  to  the  seenario  when  they  are  moving. 

We  represent  the  target  loeation  and  veloeity  as 

'P  =  [Px,Py^Pz]^  (1) 

^  —  ['^x  '^y  ^  (2) 

Therefore,  we  ean  define  the  target  state  veetor 

P  =  [Px,Py,Pz,Vx^Vy^,Vz]‘  (3) 

Under  this  set  up,  the  direet  or  referenee  path  delay  between  the 
zth  transmitter  and  the  jth  reeeiver  is 


Similarly,  the  delay  and  Doppler  eorresponding  to  the  target 
eeho  path  between  the  same  transeeiver  pair  ean  be  given  as 
[26] 

^^■  =  “(ll^“T||  +  iiV-^ll),  (5) 

Jd.,  =7  (V,  Vtj)  ,  (6) 


where  c  denotes  the  speed  of  wave  propagation  in  the  medium 
and  fc  denotes  the  earrier  frequeney;  T? ^ denote  the  unit 
veetor  from  the  zth  transmitter  to  the  target  and  the  unit  veetor 
from  the  target  to  the  jth  reeeiver,  respeetively;  ||.||  represents 
the  I2  norm;  and  ( , )  is  the  inner  produet  operator.  Clearly,  the 
dependenee  of  the  delays  and  Doppler  shifts  on  the  multistatie 
geometry  is  non-linear  in  nature. 

As  in  [15],  we  eonsider  downlink  UMTS  signals  as  our  illu¬ 
minators  of  opportunity.  Let  Ui{t)  denote  the  baseband  signal 
eorresponding  to  the  ith  transmitter.  Then 

1 

Ui{t)  =  ^  Ci„gi{t  -  nT),  (7) 

^  n=0 


where  Ci^  are  the  transmitted  quadrature  phase  shift  keying 
(QPSK)  symbols  from  the  zth  transmitter.  N  is  the  number  of 
symbols  and  T  is  the  inverse  of  the  ehip  rate.  All  the  QPSK 
symbols  are  independent  and  identieally  distributed  aeross  both 
the  transmitter  and  symbol  indiees. 


1,  i  —  i'  —  n'\ 
0,  otherwise. 


(8) 


The  pulse  gi{t)  is  defined  as  delayed  root-raised-eosine 
(RRC)  pulse  gi{t)  =  /z^  (t  —  y ) ,  where  the  delay  D  is  seleeted 
as 

o 

2 

/  h^{t)dt  =  0.95. 

-D/2 


The  RRC  pulse  h{t)  is  defined  as  follows 

hi{t) 


sin  ((^)  (1  -  aj))  +  (^)  cos  ((^)  (1  -f  gQ) 


f  ttZ  5 

i  (  Aait  \  ^ 

vM 

^  V  T  J 

h{t  =  0) 


The  roll-off  faetor  ai  lies  in  the  interval  [0, 1  ]  for  all  transmit¬ 
ters.  Further  analysis  on  the  properties  for  these  waveforms  for 
different  parametrie  values  are  give  in  [15].  We  skip  it  here  be- 
eause  it  is  not  the  main  foeus  of  our  work.  We  ehoose  ai  =  0.22, 
T  =  0.26  /zs,  and  NT  =  0.1  s.  These  are  typieal  values  for 
UMTS  systems.  Having  presented  the  transmit  signal  model  and 
the  radar  geometry,  we  now  move  on  to  the  two  different  oper¬ 
ating  seenarios  in  the  following  seetions. 


III.  Non-Coherent  Processing 
A.  Measurement  Model 

Non-eoherent  proeessing  is  useful  to  fully  exploit  the  spatial 
diversity  gain.  Under  this  operating  mode,  we  assume  the  target 
eonsists  of  several  individual  isotropie  seatterers,  similar  to 
the  model  in  [17].  We  assume  the  target  RCS  varies  with  the 
angle  of  view  and  the  different  antennas  are  sufiieiently  spaeed 
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apart  to  decorrelate  these  RCS  values.  In  other  words,  these 
RCS  values  are  uncorrelated  zero-mean  complex  Gaussian 
random  variables.  If  the  target  RCS  values  corresponding  to 
certain  bistatic  pairs  are  weak,  it  is  highly  likely  that  they  will 
be  compensated  for  by  other  pairs  with  strong  returns.  Let  the 
target  RCS  corresponding  to  the  zth  transmitter  and  the  jth 
receiver  be  denoted  by  pij .  Using  a  narrowband  assumption  on 
the  waveforms,  we  can  express  the  received  signal  at  the  jth 
receiver  after  down  conversion  from  carrier  frequency  as 


M'-p 


j27T(f^  — 


i=l 


Expressing  in  terms  of  the  RRC  pulses,  we  have 


Mt  ^  N-l 

yjW  =  i  -  nT 

i=l 


2 


j27T 


xe 


+  nij{t). 


The  reference  or  direct  path  signal  from  the  zth  transmitter  to 
the  jth  receiver  is  given  as 

yfjii)  =  Ui{t-  T^) 

nij(t)  and  nfj(t)  denote  the  additive  noise  terms.  We  assume 
that  all  the  direct  path  and  echo  signals  from  different  trans¬ 
mitters  are  separable  at  each  receiver.  This  is  a  reasonable 
assumption  as  the  different  transmitters  in  a  UMTS  system 
would  not  interfere  with  one  another  by  using  separate  spectral 
resources.  Next,  we  will  describe  the  target  detector  for  this 
problem.  We  will  use  this  detector  to  compute  the  ambiguity 
function  expression. 

B.  Matched  Filter  Based  Detector 

Given  this  measurement  model,  the  matched-filter  based  de¬ 
tector  reduces  to  computing  the  non-coherent  test  statistic  (log- 
likelihood)  corresponding  to  a  state  vector  [17],  [27] 

M'j'  Mu 

cx  (9) 

i=l  j  =  l 

where 

ytj  =  f  yjWyfj*  -  '^ij  +  ’^ti) 

Expressing  the  reference  signal  in  terms  of  the  transmitted 
signal,  we  obtain 

yi3=Jyj(^')  (y'^i  +  ?%•''(*)) 

The  noise  term 

=  nfj*  {t  -  +  rfj)  . 

Please  note  that  the  test  statistic  does  not  depend  on  the  phase 
term  and  there  are  no  interference  terms  from  the 

other  transmitters.  Even  e  does  not  influence  the 

test  statistic.  Therefore,  phase  synchronization  is  not  necessary 
between  transmitters  and  receivers  for  this  detector.  Even 
though  matched  filter  based  detector  is  the  most  well-known 


and  a  simple  way  to  approach  this  problem,  no  claims  can  be 
made  about  optimality  because  the  detection  problem  is  no 
longer  a  point  hypothesis  testing  due  to  the  noisy  nature  of  the 
reference  signal  and  instead  becomes  a  composite  hypothesis 
testing  problem.  Other  more  complex  signal  processing  ap¬ 
proaches  like  GERT  may  be  used  but  again  optimality  cannot 
be  claimed.  However,  note  than  when  the  signal  to  noise  ratio 
of  the  direct  path  or  reference  signal  is  high,  asymptotically  the 
matched  filter  detector  does  tend  towards  optimality  because 
the  reference  signal  will  tend  towards  the  transmitted  signal. 


C.  Ambiguity  Function 

We  compute  the  average  ambiguity  function  corresponding 
to  two  target  states  ji  =  [pxjPy^Pz^Vx^Vy.Vz]  and  p'  = 

\Px>Py>P'z^K^'Uy>K]  as  or  A{p^,Py,p^,V^,Vy,V^,p'^, 

Py^p'^^v'^^Vy^v'^).  fi  denotes  the  true  target  state  and  we  are 
are  comparing  with  the  state  p' .  The  averaging  is  done  over 
the  transmitted  QPSK  symbols.  As  we  mentioned  earlier,  these 
symbols  are  independent  across  the  transmitter  and  symbol 
indices.  Using  the  above  detector,  we  define  the  average  ambi¬ 
guity  function 


f  -1  Mt  Mr 

z=l  J  =  1 

where 

= 

'Ttj 

J  Ui{t-  e  (t  - 

(10) 


i27r/^,  t 


dt. 


Substituting  the  expressions  for  the  transmitted  waveforms,  we 
get 


^  ly  —  ±  ly  —  ±  « 

'  =  ^1:1:/ 


xe 


n=0  n'=0 


Cir^hi  [t  -  nT  -T^- 


-  Cin^K 


I'T 


D 


ru  I  e  dt. 


Define  as  the  monostatic  ambiguity  function  corre¬ 

sponding  to  the  individual  RRC  pulses  from  the  Ah  illumi¬ 
nator.  Using  the  expressions  we  derived  in  the  Appendix  for 
/  2 

E{\(j)^b  I  },  we  have  the  average  ambiguity  function 


— 


1 


MtMrA2 


EE 

i=l  j  =  l 


sin  (^ATTA/^f') 


+ 


sin 

AT-l 

E  (^-H) 


X,.  (Ar^.^',A/Sf;) 


n^0,n=  — N  +  l 


The  values  of  (Ar/^-^  ?  easily  computed 

using  the  result  on  monostatic  ambiguity  function  in  [15].  Note 
that  the  dependence  of  bistatic  delay  and  Doppler  terms  on  the 
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target  positions  and  velocities  is  non-linear  as  we  observe  from 
the  expressions  in  Section  11.  Further,  the  ambiguity  function  de¬ 
pends  on  the  multistatic  geometry  of  the  passive  radar  system.  In 
other  words,  an  ambiguity  function  centered  around  the  origin  is 
no  longer  sufficient  to  completely  assess  the  performance  of  our 
system.  There  are  two  parts  in  the  average  ambiguity  function 
summation  we  presented  above.  While  the  first  part  denotes  the 
self-ambiguity  of  the  pulses,  the  second  part  denotes  the  impact 
of  the  interaction  between  the  different  pulses  on  the  ambiguity 
function. 

From  [15],  we  observe  that  the  monostatic  ambiguity  func¬ 
tions  corresponding  to  the  individual  RRC  pulses  can  be  ex¬ 
pressed  as 

x,. 


oo 

=  J  (t  -  )  e 

—  OO 

=  Hi  (A/gf')  ®  H,  (A/Sf')  e^' 


dt, 


where  )  denotes  the  Fourier  transform  of  the  RRC 

pulse  hi  (t)  and  is  the  convolution  operator.  Using  the  above 
expression,  when  the  delay  =  0,  the  term  representing 

the  inter  pulse  interaction  becomes 


Therefore,  we  compute  the  zero-Doppler  cut  for  the  non-co- 
herent  ambiguity  function  [see  the  equation  at  the  bottom  of  the 
page].  The  zero-delay  cut  of  this  ambiguity  function 
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where  \Hi{Af!^^, )  (g)  )|  is  computed  in  equation  (18) 

of  [15].  Using  numerical  simulations,  we  will  plot  the  zero- 
delay  and  zero-Doppler  cuts  computed  above  to  study  the  am¬ 
biguity  performance,  later  in  the  paper. 


IV.  Coherent  Processing 
A.  Measurement  Model 

In  the  non-coherent  processing  mode  presented  above,  the 
target  induces  uncorrelated  phase  shifts  across  the  different 
bistatic  transmitter-receiver  pairs  due  to  the  widely  separated 
nature  of  the  antennas.  As  a  result,  it  is  impossible  to  coher¬ 
ently  combine  the  outputs  at  the  different  receivers.  Now, 
we  consider  a  different  scenario  under  which  such  coherent 
processing  is  possible  [16].  In  this  set  up,  the  target  consists 
of  a  single  point  scatterer  that  has  an  isotropic  complex  re¬ 
flectivity  denoted  by  (3.  This  arrangement  preserves  the  phase 
information  across  the  different  bistatic  pairs  since  there  is  no 
multiple  scatterer  interaction.  From  now  on,  when  we  mention 
to  coherent  processing,  we  refer  to  this  target  model.  Note  that 
perfect  synchronization  (phase,  time,  and  frequency)  to  ensure 
coherent  processing  is  possible  only  when  the  antennas  are 
not  too  far  apart.  Otherwise,  we  have  to  resort  to  non-coherent 
processing  described  in  the  previous  section.  In  this  paper,  we 
are  computing  the  expressions  for  both  these  cases. 

Under  a  narrowband  assumption,  the  baseband  received 
signal  at  the  jth  receiver 
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The  reference  path  signal  yj  remains  the  same  as  in  the  non¬ 
coherent  processing  mode  since  the  only  difference  between  the 
two  scenarios  stems  from  the  target  echo 
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B.  Matched  Filter  Based  Detector 

The  coherent  test  statistic  for  a  matched  filter  based  detector 
corresponding  to  the  target  state  /i  is  given  as 


L^{y)  oc 


i=l  j=l 


(11) 


where 

yij  =  /  +  '^fj) 


We  substitute  the  expression  for  the  reference  signal  to  get 

ytj  =  / yj(*^ 

-j2Tr( 

xe 

where 


(0  =  nfj*  {t  -  t[j  +  T^)  e 


The  dependence  of  this  detector  on  the  phase  term  is  evident 
from  the  test  statistic. 

C  Ambiguity  Function 

For  this  case,  the  ambiguity  function  must  reflect  the  coherent 
combination  of  the  different  bistatic  pair  components  instead  of 
just  combining  the  amplitudes  or  energies.  Therefore,  we  define 
the  average  coherent  ambiguity  function  as  [17],  [27] 
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After  expanding  the  absolute  value  of  the  coherent  summation, 
we  can  express  the  average  coherent  ambiguity  function  as 
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Using  the  derivations  we  made  in  the  Appendix,  we  compute 
the  coherent  ambiguity  function 
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Using  the  results  in  [15],  it  can  be  easily  verified  that  the 
zero-delay  cut  of  the  coherent  ambiguity  function 
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Fig.  2.  \xh(nT  +  t,  0) |  for  different  n  as  a  function  of  r. 


Fig.  3.  \xh{nT,  u) \  for  different  n  as  a  function  of  i/. 


Similarly,  the  zero-Doppler  cut  of  the  coherent-ambiguity 
function 
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Note  that  unlike  the  non-coherent  scenario,  it  was  not 
possible  to  express  the  average  ambiguity  function  for  the 
coherent  case  as  a  summation  of  just  the  amplitude-squared 

terms  ,  A/^^,  )|  due  to  the  presence  of  the  phase 

terms.  Therefore,  we  expressed  in  terms  of  the  complex  mono¬ 
static  ambiguity  function.  The  second  summation  term  in  the 
average  ambiguity  expressions  above  represents  the  interaction 
in  between  the  different  pulses  of  different  symbol  intervals. 
Further,  even  in  the  coherent  case,  the  ambiguity  function 
depends  on  the  multistatic  radar  geometry,  i.e.,  the  locations  of 
the  transmitters,  receivers,and  the  target  cell  of  interest. 

V.  Numerical  Examples 

We  first  examine  the  impact  of  the  terms  denoting  the  inter¬ 
action  between  the  pulses  in  the  zero-Doppler  and  zero-delay 
cuts,  i.e.. 
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for  varying  n.  From  Fig.  2,  we  observe  that  the  cross-terms  do 
have  peaks  around  r  =  —nT.  This  is  because  at  those  de¬ 
lays,  the  two  pulses  in  the  ambiguity  function  from  different 
symbol  exactly  overlap.  However,  the  scaling  factor  will 

attenuate  these  peaks  heavily.  Since  the  integration  time  for 
passive  radar  is  typically  large,  the  value  of  N  causes  severe 
attenuation  to  these  cross-pulse  peaks  making  them  insignifi¬ 
cant  in  computing  the  ambiguity  function.  From  Fig.  3,  we  no¬ 
tice  that  the  impact  of  the  cross-pulse  terms  is  more  significant 
at  very  high  frequencies.  However,  typical  Doppler  frequen¬ 
cies  are  present  at  the  lower  end  of  the  spectrum.  Around  the 
zero-Doppler  point,  the  contribution  of  these  cross-pulse  terms 
is  negligible.  This  contribution  is  further  attenuated  when  we 
multiply  by  .  Therefore,  we  ignore  these  terms  for  the 

simulations. 
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Fig.  5 .  Zero-Doppler  cut  of  non-coherent  ambiguity  function  when  [p^ ,  Py  ]  = 
[3,  4]  km  and  [v^,  t’y][30,  50]  m/s. 


We  use  numerieal  simulations  to  demonstrate  the  ambiguity 
performanee  of  distributed  passive  multistatie  radar  systems. 
For  simplieity,  we  eonsider  the  targets  and  antennas  to  be 
present  in  the  x-y  plane  but  the  results  will  hold  true  even  for 
the  three  dimensional  ease.  We  eonsider  5  transmitters  and  5 
reeeivers  (See  Fig.  4)  loeated  at 


ti  =  [0,  2]  km; 
^2  =  [0,  3]  km; 

=  [0,  5]  km; 
^4  =  [2,  2]  km; 
=  [6,0]  km; 


Ti  =  [3,  0]  km; 

=  [5,  0]  km; 
rj  =  [0,3.5]  km; 
rt  =  [1, 1]  km; 

=  [4,  4]  km. 


We  eompute  the  ambiguity  funetion  in  the  Cartesian  domain 
eentered  around  the  position  [3,4]  km  and  veloeity  [30,50]  m/s. 
We  ehose  the  same  system  parameters  as  in  [15]  for  the  simu¬ 
lations;  observation  time  NT  =  0.1  s,  T  =  0.26  /xs,  a  =  0.22, 
and  the  eenter  frequeney  fc  =  2100  MHz. 


Fig.  6.  Non-coherent  ambiguity  profile  in  both  position  dimensions. 
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Fig.  7.  Zero-delay  cut  of  non-coherent  ambiguity  function  when  [px,Py]  = 
[3,  4]  km  and  [v^,  Uy][30,  50]  m/s. 


A.  Non-Coherent  Processing 

First,  we  begin  with  the  non-eoherent  operating  seenario.  In 
Fig.  5,  we  plot  the  zero-Doppler  eut  of  the  ambiguity  funetion. 
We  elearly  observe  that  the  peak  eorresponds  to  the  true  target 
position  that  we  are  testing.  Note  that  unlike  a  monostatie  or 
a  bistatie  eonfiguration,  we  get  a  unique  peak  for  the  two  di¬ 
mensional  ambiguity  funetion  instead  of  a  peak  eirele  or  ellipse, 
respeetively.  Further,  when  multiple  targets  are  present,  mul¬ 
tiple  peaks  eorresponding  to  different  targets  will  be  present  at 
the  final  output  of  the  matehed  filter  proeessing.  In  order  to  ob¬ 
serve  the  resolution  in  eaeh  dimension,  we  further  fix  the  target 
loeation  in  eaeh  dimension  to  the  true  value  and  plot  the  am¬ 
biguity  profile  in  the  other  dimension  (see  Fig.  6).  We  define 
resolution  as  the  distanee  between  the  20%  drop-off  points  in 
the  main  lobe.  We  observe  that  the  level  drops  by  20%  in  the 
x-dimension  at  2850  m  and  3150  m.  Therefore  the  position  res¬ 
olution  in  the  x  dimension  is  approximately  300  m.  However, 
in  the  y-dimension,  the  level  drops  below  20%  at  3930  m  and 
4070  m,  respeetively.  Therefore,  the  resolution  is  mueh  better  in 
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Fig.  8.  Non-coherent  ambiguity  profile  in  both  velocity  dimensions. 


Coherent  Ambiguity  Function  A 


1 


Fig.  9.  Zero-Doppler  cut  of  coherent  ambiguity  function  when  [px,Py]  = 
[3,  4]  km  and  [v^ ,  [30,  50]  m/s. 

the  y  dimension  (approximately  140  m).  Note  that  these  resolu¬ 
tions  depends  on  the  geometry  of  this  speeifie  example  and  will 
ehange  depending  the  loeations  of  the  transmitters,  reeeivers, 
and  targets. 

Next,  we  plot  the  zero-delay  eut  of  the  ambiguity  funetion  in 
Fig.  7.  The  peak  eorresponds  to  the  true  values  of  the  target  ve- 
loeities  that  we  are  testing,  i.e.,  [30,50]  m/s.  Just  as  we  did  for 
the  target  positions,  we  fix  the  target  veloeity  in  eaeh  dimen¬ 
sion  and  eompute  the  ambiguity  profile  in  the  other  dimension 
(see  Fig.  8).  The  veloeity  resolutions  in  the  x  and  y  dimensions 
are  approximately  5.5  m/s  and  2.5  m/s,  respeetively.  Therefore, 
both  the  position  and  veloeity  resolutions  are  better  in  the  y  di¬ 
mension  for  this  example. 

B.  Coherent  Processing 

Now,  we  present  the  eoherent-ambiguity  funetion  for  the 
same  example  as  the  non-eoherent  ease.  We  expeet  the  ad¬ 
ditional  information  in  the  form  of  phase  will  give  enhaneed 
resolution  performanee  when  eompared  with  non-eoherent 
proeessing.  First,  in  Fig.  9,  we  plot  the  eoherent  zero-Doppler 


Fig.  10.  Coherent  and  non-coherent  ambiguity  functions  in  the  x-position 
dimension. 


Py  (in  m) 

Fig.  11.  Coherent  and  non-coherent  ambiguity  functions  in  the  y-position 
dimension. 

eut.  We  notiee  that  the  ambiguity  funetion  is  mueh  sharper 
when  eompared  with  the  non-eoherent  ease.  Please  note  that 
we  have  zoomed  into  the  ambiguity  plot  by  ehanging  the  seales 
on  the  X  and  y  axes  to  eapture  the  exaet  resolution. 

We  fix  the  y-position  to  the  true  value  and  plot  the  ambiguity 
funetion  along  the  x-dimension  in  Fig.  10.  We  observe  that  the 
resolution  for  the  eoherent  ease  is  0.06  m  as  opposed  to  300  m 
for  non-eoherent  ease.  We  used  the  same  20%  drop-off  defini¬ 
tion  for  the  resolution.  In  Fig.  11,  we  plot  the  ambiguity  fune- 
tions  in  the  y-dimension.  The  eoherent  resolution  is  0.1  m  eom¬ 
pared  to  140  m  using  non-eoherent  proeessing.  Even  though  the 
resolution  is  slightly  poorer  in  the  y-dimension  when  eompared 
with  the  x-dimension,  we  observe  that  the  side-lobes  immedi¬ 
ately  surrounding  the  main  lobe  are  suppressed  to  a  larger  extent 
in  the  y-dimension  for  this  example.  Note  that  even  though  there 
are  signifieant  sidelobes  in  the  eoherent  proeessing  mode,  the 
peaks  of  these  sidelobes  still  have  a  mueh  lower  value  than  the 
eorresponding  ambiguity  on  the  mainlobe  of  the  non-eoherent 
proeessing  mode. 
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Coherent  Ambiguity  Function  ^ 


(in  m/s)  so 


Fig.  12.  Zero-delay  cut  of  coherent  ambiguity  function  when  [px,Py]  = 
[3,  4]  km  and  [v^,  Uy][30,  50]  m/s. 


VI.  Concluding  Remarks 

We  considered  a  passive  distributed  radar  system  using 
the  UMTS  downlink  as  the  signals  of  opportunity.  First,  we 
described  the  multistatic  measurement  model  for  this  system. 
Next,  we  presented  the  matched  filter  based  detectors  under 
two  different  operating  scenarios;  non-coherent  and  coherent. 
Then,  we  computed  the  ambiguity  functions  under  both  these 
scenarios  to  analyze  the  performance.  We  used  numerical 
examples  to  demonstrate  our  analytical  results.  We  observed 
that  passive  multistatic  radar  systems  offer  sharp  resolutions 
when  operating  in  coherent  processing  mode. 

Since  the  ambiguity  analysis  we  derived  in  this  paper  is 
geometry-dependent,  in  future  work,  we  will  use  these  results 
along  with  our  Cramer-Rao  lower  bound  analysis  in  [31]  to 
optimally  select  the  illuminators  of  opportunity  from  amongst 
several  possible  choices.  Further,  we  will  extend  our  analysis 
beyond  UMTS  signals  to  include  other  classes  of  illuminators. 

Appendix 
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Fig.  13.  Coherent  ambiguity  profile  in  both  velocity  dimensions. 

Fig.  12  plots  the  zero-delay  cut  of  the  coherent  ambiguity 
function.  We  observe  that  the  peak  coincides  with  the  true  target 
velocity  values.  In  Fig.  13,  we  plot  the  coherent  ambiguity  func¬ 
tion  in  both  the  velocity  dimensions  by  fixing  all  the  other  pa¬ 
rameters  to  the  true  values.  We  noticed  that  the  resolutions  for 
the  coherent  case  are  approximately  1 .5  m/s  and  1 .6  m/s  in  the  x 
and  y  dimensions,  respectively.  Clearly,  the  resolution  has  im¬ 
proved  when  compared  with  the  non-coherent  case. 

Note  that  the  coherent  and  non-coherent  ambiguity  functions 
represent  the  matched  filter  based  detectors  under  two  different 
target  models.  Therefore,  even  though  coherent  processing 
offers  the  enhanced  resolution,  it  is  not  applicable  when  the 
target  consists  of  multiple  individual  isotropic  scatterers, 
leading  to  uncorrelated  phase  shifts  across  different  bistatic 
pairs.  Further,  phase  synchronization  between  all  the  transmit¬ 
ters  and  receivers  is  critical  for  coherent  processing  and  this 
may  not  be  possible  always  due  to  several  physical  limitations 
like  inaccurate  knowledge  of  the  antenna  locations  and  local 
oscillator  characteristics  [28]-[30]. 
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For  the  non-coherent  processing  scenario,  we  need  to  com¬ 
pute  the  amplitude-squared  terms 
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Since  the  QPSK  symbols  are  independent  at  different  symbol 
periods, 

1,  n  —  n’  ^  k  — 

E  {^inCin'Cik^^ik']  =  ^  1,  n  =  k,  n'  ^  n,  n'  =  k'; 

0,  otherwise. 


Let  denote  the  monostatic  ambiguity  function  corre¬ 

sponding  to  the  individual  RRC  pulses  from  the  ith  illuminator. 
The  summation  can  be  split  into  these  two  conditions  as  below 
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§  s  /  / 
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.  .i' 
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N-1  N-1 

E  E 

n=0  k' ^n,k'  =0 


X  e  ^^3  f 


hdt  -  nT  -  —  - 


^  fl'  \ 

y  -  j  e 


xhi  (t'  -k'T-j-  tA  I  . 


■g^^'  -j27r  I  /c  Arf/  +/£  r^.  -f^  r^. 

g  j  =  6  j 


J2^  (/cAry; 


Expressing  the  transmitted  waveforms  in  terms  of  the  indi¬ 
vidual  eonstituent  RRC  pulses,  we  obtain  the  following 


E{€r'*€^} 


N-1  N-1  N-1  N-1 


^{^inC^n'Ci^k4k'} 


n=0  n'=0  k=0  k'=0 


A/r.  =fL-fL 


r  ,  a'l  sin  ( nNT A fZ)^ 

}=  -(rAr') 

J  sm  TrTA/gM 


2  iv-1  iv-i  2 

A/Dr')|  +E  E  |x/».(Ar^'''+(*'-n)r,A/gf')  . 

n=0  k' ^n,k'  =0 

Note  that  the  seeond  part  of  the  summation  denotes  the  impaet 
of  the  interaetion  between  the  different  pulses  on  the  ambiguity 
fimetion. 


X  hi  n'T  “  y  “  ^  * 

xy, 

X  4  U  -  k'T  -  —  -  4-,  j  e  dtdtk 


Sinee  the  transmitted  QPSK  pulses  are  independent  aeross  the 
different  transmitter  indiees  and  symbol  intervals,  we  have 

'1,  n  =  k  =  k'; 

T-ir*  *  d  _  I5  TL  =  k^  n  ^  71, 

^  i^in^in' Ci' kCi' k' }  —  ^  . 

<  0,  otherwise. 


N‘^E  I 


sin  (ttA^TA/^^,  )  ,  ,  ,.2 

sm(7rTA/SfJ  ^ 

+  ^  {N-\n\)\xh.{ATl^'^'  +  nT,Af^^)\\ 

n^0,n=  —  N-\-l 

For  the  eoherent  proeessing  ease,  instead  of  the  amplitude- 
squared  terms,  we  need  to  eompute  the  following  terms 


=  e  j 
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X  e  ^3 
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xe  u*,it'-Tp,.,)e  At' 


Substituting  the  above  expression,  we  obtain 
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N‘^E 
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%  —  i\ 

0,  otherwise. 

n=0  k=0 
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+  6u'(j:  E 
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where  [see  the  equation  at  the  top  of  the  page].  Upon 
simplifieation, 
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Substituting  this  expression  gives  the  average  ambiguity  fune- 
tion 
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